Andreev Bound states as a phase sensitive probe of the pairing symmetry of the iron 

pnictide superconductors. 
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A leading contender for the pairing symmetry in the Fe-pnictide high temperature superconductors 
is extended s-wave s±, a nodeless state in which the pairing changes sign between Fermi surfaces. 
Verifying such a pairing symmetry requires a special phase sensitive probe that is also momentum 
selective. We show that the sign structure of s± pairing can lead to surface Andreev bound states 
at the sample edge. In the clean limit they only occur when the edge is along the nearest neighbor 
Fe-Fe bond, but not for a diagonal edge or a surface orthogonal to the c-axis. In contrast to d- 
wave Andreev bound states, they are not at zero energy and, in general, do not produce a zero 
bias tunneling peak. Consequences for tunneling measurements are derived, within a simplified two 
band model and also for a more realistic five band model. 

PACS numbers: 



Central to the physics of the recently discovered iron 
pnictide superconductors, is the nature of the pairing 
symmetry. Given the high transition temperatures and 
proximity to a magnetically ordered state, a non-phonon 
mechanism is indicated. Hence, unconventional pairing 
symmetry, as often occurs in other such cases like the 
cuprates, organic conductors and heavy fermion materi- 
als, is to be expected. Current theories predict a wide 
variety of possible pairing states, ranging from p-wave^ 
to d-wave&i, and also extended s-wave s±, a nodeless 
pairing state where the pairing changes sign between dif- 
ferent Fermi pockets^&i^. The last is a theoretically 
attractive candidate since it is a naturally connected 
with the magnetic order observed in the undoped par- 
ent materials&£&^. While initial experimental results 
seemed to favor a pairing state with nodes in the gap 
function, recently, the bulk of the experimental evidence 
has been in favor of a completely gapped state. For ex- 
ample, ARPES experiments— see a full gap around all 
Fermi pockets, and recent specific heat measurements are 
consistent with the absence of nodesii. The extended 
s-wave state is therefore emerging as an strong candi- 
date for the pairing symmetry. Hence, it is important 
to consider experimental tests that can directly probe 
such a state. In contrast to the d-wave symmetry of the 
cuprates, which changes sign under rotation and hence 
can be probed by a class of phase sensitive experiments^, 
the sign change in the s± state is not accessible by a 
real space transformation. Instead, as shown in figure 
l(a)[ the pairing phase changes on moving across the 
Brillouin zone. Probing such a sign change is a much 
more challenging task. A parametrization of the pair- 
ing function that gives this sign structure on the fermi 
surfaces is A(k x ,k y ) = Aq cos k x cos k^, where (k x ,k y ) 
represents the crystal momentum along the Fe-Fe bonds. 
The pairing changes sign on crossing the k x = ±7r/2 and 
k y = ±7r/2 lines. It is this sign change that we wish to 
probe. 

An interesting consequence of d-wave pairing in the 



cuprates is the presence of Andreev bound states at 
the edge which leads to a zero-bias peak in conductiv- 
ity in superconductor/normal junction o 13 i 14 ' 15 . The ex- 
istence of these modes can be intuitively understood as 
follows. Consider an electron with momentum k inci- 
dent on a surface, where it is reflected to momentum 
k' . When the pairing function changes sign in this pro- 
cess A(k) / A(k') < 0, then the problem is related to a 
superconductor-normal-superconductor junction with a 
relative phase shift of ir. There it is well known that a 
midgap state arises near the chemical potential^ 6 . This 
leads to Andreev bound states along those edges of the 
d-wave superconductor that reflect particles in the sign 
reversed way described above. Note, the Andreev bound 
states reach all the way down to the chemical potential 
(zero energy), and hence are experimentally detected as 
a zero bias peak in tunneling measurements^^. The 
dependence on edge orientation can be erased in sam- 
ples with disordered surfaces. Nevertheless, in some 
cases a strong experimental difference has been observed 
between 45° rotated a-b plane edges^. The zero bias 
anomaly does not arise in c-axis tunneling, as predicted. 

In this letter we study the surface Andreev bound 
states at the edge of an extended s-wave superconductor. 
Although our results are more general, we focus on the 
models proposed for the pnictide superconductors. By 
the intuitive argument described above, we expect that 
when an electron is scattered by the edge between Fermi 
pockets with opposite sign of pairing, Andreev bound 
states should result. We indeed find this to be the case, 
in the microscopic models that we study, but with one 
important distinction from the d-wave case. The spec- 
trum of Andreev bound states within the bulk gap does 
not reach down to zero energy. A zero bias tunneling 
anomaly is hence not expected. Nevertheless, the pres- 
ence of these Andreev modes will lead to peaks in the 
tunneling current within the bulk gap, which should be 
observable. From figure 1(a) it is seen that scattering 



between Fermi surfaces with opposite sign of pairing oc- 
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curs for an edge parallel to the nearest neighbor Fe-Fe 
bond (along the [100] or [010] directions) but not for the 
diagonal edge ([110] direction) or the c-axis edge ([001] 
direction - we assume negligible c-axis dispersion). In- 
deed, we only find bound Andreev states for the expected 
edge directions, and these disappear if the relative sign of 
pairing is removed (i.e. with regular s- wave pairing) . The 
sensitivity to edge orientation is a key signature of these 
Andreev surface states. At present, tunneling experi- 
ments on the Fe-pnictide superconductors are restricted 
to polycrystalline samples^, which makes comparison to 
theory difficult. Tunneling experiments on single crystals 
can unambiguously identify these states. 



We study two microscopic models of the Fe pnictides. 
The first is a simplified two band modeli&is', in which 
the Andreev bound states can be studied analytically. 
The Andreev bound state energy if found to be con- 
trolled by the overlap between the Bloch wavefunctions 
of the incident and scattered electronic states, and sat- 
isfies the intuitive requirement that as the wavefunction 
overlap decreases, the in-gap state approaches the con- 
tinuum. With increasing overlap, the bound state ap- 
proaches zero energy, which is achieved if the overlap is 
perfect. A typical choice of parameters yields a dispers- 
ing in-gap mode whose minimum energy is less than 50% 
of the bulk gap. Since the two band model differs in some 
subtle but important respects from the LDA determined 
band structures^, we have also studied a more realistic 
five-band tight binding models. Here we find qualita- 
tively similar results. For the particular doping that we 
study in detail, the in-gap state is at a higher energy, 
92% of the bulk gap at that edge momentum. Since a 
perfectly smooth edge is hard to realize experimentally, 
we model edge roughness by allowing for orbital mix- 
ing during scattering (but preserve the edge momentum 
quantum number). Within this crude model of a rough 
edge, the in-gap states occur at much lower energy (eg. 
50% of the bulk gap value). Observable signatures in 
experiments with planar and point contact tunnel junc- 
tions are discussed. In the treatment below, we solve 
the Bogoliubov-deGennes equations for superconductor 
quasiparticles near the sample edge, with a uniform pair- 
ing potential. For weak coupling superconductors, the 
pair potential is to be calculated self consistently, which 
we do not do here. A self consistent solution is expected 
to lead to a smaller pairing potential near the edge, due 
to the presence of in-gap states. This is not expected to 
change our results qualitatively, and moreover, it is not 
even clear if such weak coupling reasoning is justified for 
these high temperature superconductors. 



Two Band Model: We first present the the study of 
the two band model for pnictides at half filling. The 
details of this model could be found elsewhere^. The 
band Hamiltonian we consider in momentum space has 
the form: 
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FIG. 1: a) Fermi surfaces of Fe pnictides and extended s- 
wave pairing sign in the unfolded Brillouin zone b) Real space 
square lattice of iron atoms and the [100] boundary (red) 
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where ^1(^2) represents 3d xz (3d yz ) orbitals of Iron, e 
is the 2x2 antisymmetric matrix and a is the spin index. 
The momentum sum is over the enlarged Brillouin zone 
k x E [— 7T, n), k y 6 [— 7r, 7r). The 2x2 matrix K(k) in 
Ref^a 



after rotation into the d xz ,y Z basis is: 
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and "71,2,3 are Pauli matrices and / is the identity matrix. 
Here we present the numerical results for t\ = leV ,t 2 = 
1.7 'eV and t{ = t' 2 = at zero doping^. 

We numerically diagonalize a 600 site wide lattice, with 
an edge along the [100] direction. The momentum fcy 
along this direction is conserved and used to label the 
eigenstates. Interestingly we see (Figure [U that a dis- 
persing mid-gap state is obtained, with energy as low 
as 39% of the bulk gap. This is the Andreev surface 
bound state. For convenience, we work with a bulk gap 
of 46 meV which is a few times larger the experimentally 
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measured gap. We have checked that the energy of the 
mid-gap features scale linearly with the gap size over a 
wide range hence our results hold when scaled with the 
bulk gap. The results were also compared against those 
obtained from transfer matrix method 2 -! which does not 
suffer from finite size errors. The dispersing mode leads 
to a feature in the density of states (Figure inset), 
which could potentially be detected in the tunneling ex- 
periments. No surface Andreev states are obtained with 
the edge along the [110] direction, or if an onsite s-wave 
pairing, with uniform sign in momentum space is used. 
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FIG. 2: Dispersion of Andreev bound states (red) and the 
bottom of the continuum (black) as a function of feii , the mo- 
mentum parallel to the edge. Energy is scaled by A, the 
minimum energy gap in the bulk. Inset shows the density of 
states integrated over feii. 

In the two band model, one feature of the edge mode 
as can be seen in figure [2] is that at transverse momen- 
tum fen = the edge mode merges with the continuum. 
This arises because the Bloch wavefunctions on the elec- 
tron and hole pockets at fen =0 are orthogonal (they 
have opposite eigenvalues under reflection). Hence, the 
gap sign changing scattering is forbidden, and no sur- 
face Andreev state appears. On moving away from this 
special momentum, the scattering is allowed, leading to 
Andreev bound state. Our analytical treatment below 
captures the relation between wavefunction overlap and 
Andreev bound state energy. 

We now analytically demonstrate the presence of edge 
modes as a result of change of sign of order parameter on 
the two Fermi surfaces in the two band model. For each 
value of fen, the momentum parallel to the edge, we have 
a one dimensional Hamiltonian. The low energy excita- 
tions occur close to the Fermi points, labeled as left mover 
and right movers (fig. [3]) for each Fermi pocket. The low 
energy electronic field *& a (x) written as a two component 
object to account for the two orbitals, can be expanded 
as ¥„(x) = E„</Wz)u« e* k f- x +i> Lna (x)ut e - lk f^, 
where n labels different Fermi pockets, a labels physical 
spin (we suppress k\\ index from now on for simplicity) 



and u n are the two component Bloch wave functions 
at the Fermi points. 

R R L L 



FIG. 3: Effective one dimensional model of the two band 
Hamiltonian, at a fixed ku. Low energy excitations are la- 
beled as right (R) and left (L) movers. Hole (electron) Fermi 
pockets are shown in green (red). 

In these variables, the superconducting Hamiltonian is: 




Here n = 1, 2 are the hole and electron pockets, while 
v n labels the fermi velocity and A„ is the superconduct- 
ing gap on the n'th fermi pocket. 

We are looking for a localized solution on the edge and 
consider the following ansatz to represent the edge mode 
which diagonalize Hamiltonian [5] 

TO •* 

+ a nLl^nLl( X ) U n*) 

*2 = E / dx e-^- ik ^(ct nR ^ nRi {x)u** 

71 J 

+a„Lf0nLT( 2; ) U n) 

They satisfy [5>, H] = E $ with E < A which leads to 
BdG equations for a's. At the boundary the particle and 
hole amplitude should separately vanish 22 . This leads to 
the following boundary equation: 

E C l" U n a «flT + C2 n U^a„L T = (6) 

TO 

J2c 2n u^*a* nRi +c ln u^*a* nLl = (7) 

TO 

These are set of equations for c's. The condition 
for existence of non-zero solution for c determines the 
bound state energy. The Bloch function on the two fermi 
surfaces can then be chosen to have the simple form: 
uf = (ri,si) T , uf = (-ri,si) T , uf = (r 2 ,s 2 ) T and 
u 2 = (— r 2 ,s 2 ) T . We can solve the BdG equations for 
the bound state in general, regardless of magnitude of 
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gap on the fermi surfaces but the expressions are greatly 
simplified if we assume that |Ai| = |A 2 | = A. We present 
the result for this case here. When AiA 2 < we get a 
pair of bound states related by time reversal symmetry 
with energy: 



E = ±A 



S2 r l + r 2 s l 



S2n - r 2 si 



(8) 



Notice the energy of these states is related to the or- 
bital structure on two fermi surfaces. With real wave- 
functions, these states occur within the bulk gap if 
riS\T2S2 < 0, which is satisfied in the two band model 
of Eq. [H Only with complete orbital overlap of and 
u 2 are they at zero energy. On the other hand, when 
the gaps have the same sign (i.e. Ai = A 2 = A), no 
localized mode is found as expected. 

Five band model: We also preformed a numerical study 
of a more realistic five band model 20 at 18% electron 
doping. We again consider a Hamiltonian of the form ([2]), 
but dk,CT are now five component spinors. The detailed 
form of K(k) is given in ref. I20I . Superconducting pairing 
is as in Eqn. [31 but now with five orbitals. We performed 
exact diagonalization on lattices 600 sites wide with a 
[100] edge. While we do not present a general analytical 
solution to the five band model edge modes due to its 
complexity, we have verified that at feii = 0, where the 
problem simplifies, the mid-gap state is obtained. 
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FIG. 4: Andreev bound state dispersion (red) and bottom of 
continuum (black). Inset presents the density of states for 
transverse momenta near k\\ = (green region). Dotted lines 
denote the positions of fermi surfaces 

The dispersion of edge mode as a function of trans- 
verse momentum is shown in Fig. |4j Although Andreev 
bound states are present, in this case they happen to be 



close to the continuum. Around k 



the Andreev 



bound state is at 92% of the bottom of the continuum 
for that momentum. Note, the precise value of the edge 
mode energies depend sensitively on the band structure 
details, as well as on doping. Hence, the main value of 
this result is in establishing the existence of the Andreev 



bound states in the five band model as well. Moreover 
they also occur for fen = 0, since the wavefunction overlap 
is nonvanishing at that point, in contrast to the two band 
model. In the inset of Fig. 0] the density of states near 
fcii = is shown in units of [eV] -1 per unit-cell. This 
may in principle be probed by planar junction tunneling, 
which probes excitations in a small momentum cone. 
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FIG. 5: Dispersion of in-gap mode (red) and bottom of contin- 
uum (black) after adding a perturbing potential on the edge. 
Inset present the density of states. Dotted lines denotes the 
position of Fermi surfaces 

Up to here, we have only considered a clean edge, but 
this is not always the case in a realistic experimental set- 
ting. From our analytical treatment it is clear that the 
energy of the edge mode depends on the orbital struc- 
ture on the two fermi surfaces and so indeed any (non- 
magnetic) perturbing potential that mixes the orbitals 
on the edge can change this mid-gap energy. In Fig. [5] 
we present one example of the dispersion of mid-gap 
states when we have an impurity potential on the edge 
which transforms the orbital structure of one fermi sur- 
face to the other. Guided by this expectation we choose 
the following potential on the edge. We choose an im- 
purity potential on the boundary (x — 0) of the form 
Hi m p = Q.3S X fi[u^ <E> u 2 + h.c] where and u 2 de- 
note the Bloch wave functions for the electrons on the 
small hole Fermi surface, and the electron Fermi surface 
respectively at fen = 0, in a gauge where they are purely 
real. The bound states now reach down to about 50% of 
the bulk gap as shown in Figure^] For technical reasons, 
we do not consider the most general impurity potential 
that would break translation symmetry along the edge 
and hence we still retain the momentum fey . 

Conclusions: In this letter we have shown the presence 
of Andreev bound state on the surface of a supercon- 
ductor with extended s-wave pairing, which are strongly 
sensitive to surface orientation. In contrast to d-wave 
Andreev bound states, they are typically at a finite en- 
ergy. Although the bound state energy is hard to predict 
reliably, since it depends in detail on the orbital structure 
of Bloch wave functions on the fermi surfaces with anti- 
phase paring, the states were found to occur in different 
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microscopic models of the Fe-pnictides with extended s- 
wave pairing. Tunneling experiments on a-b plane edges 
with planar or point contacts, or even STM on step edges 
can potentially detect these states which, if observed, will 



be an important step to establishing the pairing sym- 
metry in these materials. We thank Dung-Hai Lee for 
instructive discussions. We acknowledge support from 
LBNL DOE-504108. 
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